Forcing and Velocity Correlations in a Vibrated Granular Monolayer 
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The role of forcing on the dynamics of a vertically shaken granular monolayer is investigated. 
Using a flat plate, surprising negative velocity correlations are measured. A mechanism for this 
anti-correlation is proposed with support from both experimental results and molecular dynamics 
simulations. Using a rough plate, velocity correlations are positive, and the velocity distribution 
evolves from a gaussian at very low densities to a broader distribution at high densities. These 
results are interpreted as a balance between stochastic forcing, interparticle collisions, and friction 
with the plate. 
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Granular gases, systems of large numbers of macro- 
scopic grains in rapid motion and interacting through 
dissipative collisions, appear in a wide range of indus- 
trial applications and natural phenomena. Energy must 
be supplied externally to compensate for the inelastic col- 
lisions, so granular gases are necessarily systems out of 
equilibrium. As a result, they may display dramatic non- 
equilibrium effects such as non-gaussian velocity distribu- 
tions p], 0, 3; 0, and long-range spatial velocity corre- 
lations ffl, 3, p| |10fl. Non-gaussian velocity distributions 
are a direct demonstration of the inapplicability of the 
Gibbs distribution, and significant velocity correlations 
indicate the absence of 'molecular chaos', which is a cru- 
cial approximation normally used to solve the Boltzmann 
equation and to calculate other fundamental quantities 
in kinetic theory. Recent theoretical work has focused on 
the non-equilibrium steady state obtained when the en- 
ergy supplied by spatially homogeneous random external 
forcing is balanced by the dissipation due to the collisions 
such that the average energy of the system remains con- 
stant j|, [?], ||. Non-gaussian velocity distributions and 
algebraically decaying velocity correlations arise as a di- 
rect consequence of the energy injection. Non-gaussian 
velocity distributions and velocity correlations have been 
observed in a number of experiments jj], ||, [| |l0| , but in 
each case the forcing is sufficiently different from that of 
the theoretical models that a direct comparison is diffi- 
cult. In this Letter, we provide a direct demonstration of 
the determining role that the forcing plays on the spatial 
velocity correlations in a homogeneously forced granular 
gas. We will also show how the framework of kinetic the- 
ory can provide a coherent description of the origin of 
the observed velocity distributions and correlations. 

We have investigated a quasi-2D granular system con- 
sisting of a layer of a large number of spherical parti- 
cles partially covering a vertically driven horizontal plate. 
Two plates with different surface properties were used in 
the experiments: a smooth circular plate (20 cm in di- 
ameter) made of black anodized aluminum, and a rough 
hexagonally shaped plate (30 cm between opposite cor- 
ners) . The roughness of the latter is provided by a close- 
packed lattice of blackened steel balls (1.19 mm diameter) 
glued to a flat plate. The granular gas is made of uni- 



form stainless steel spheres with diameter a — 1.59 mm 
(smooth plate), or a — 3.97 mm (rough plate). On the 
smooth plate, the layer is constrained from above by an 
anti-static coated plexiglas lid, 1.7 a above the plate. 
Using an electromagnetic shaker, the plate is driven si- 
nusoidally, and the granular layer is brought to a non- 
equilibrium steady state. The strength of the shaking 
T = Aui 2 j 1 g is directly measured with a fast response ac- 
celerometer (A is the amplitude and lj the frequency of 
the plate oscillation, and g the acceleration due to grav- 

ity). 

A high resolution camera (Pulnix TM1040) placed 
above the plate records the bright spots on the tops of 
the balls produced by a stroboscopic LED array. Instan- 
taneous horizontal velocities are obtained from the 
displacements between strobe pulses. The strobe is syn- 
chronized with the camera so that the first pulse occurs 
at the end of the exposure of one frame, and the next 
pulse occurs at the beginning of the subsequent frame. 
The time interval between pulses (~ 1 ms) is chosen to 
be significantly smaller than the mean collision time. The 
width of both pulses is typically less than 0.05 ms. Dis- 
placements of the centers of the bright spots are measured 
with 0.1 pixel accuracy. To eliminate systematic errors in 
the velocity correlations, snapshots are taken at a fixed 
phase relative to the driving signal corresponding to the 
maximum average position of the layer. 

The general behavior of the granular layer on the 
smooth plate has been reported previously m and a rich 
phase diagram depending on both T and uj was found. 
Here we focus mainly on the fluidized regime, obtained 
for r > 1. All results presented in this Letter were ob- 
tained with r = 1.5 and lo/2-k — 60 Hz. 

The longitudinal and transverse velocity correlations, 
C|| and Cj_, respectively, are calculated by 
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where the sum runs over the N r pairs of particles sepa- 



rated by a distance : 
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is the projection of along the 



line connecting the centers of particles i and j, and vf- is 
the projection perpendicular to that line. 
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FIG. 1: Longitudinal velocity correlations normalized by the 
granular temperature T on the smooth plate: F — 1.5, 
p = 0.4 (•) and p = 0.5 (o). Inset: transverse velocity cor- 
relations. Data represent averages over approximately 8000 
pairs of images with about 180 balls per image. 



Figure [l] shows Cm (inset: C±) normalized by the gran- 
ular temperature T = (v?) obtained on the smooth plate 
for p = 0.4 and 0.5, and V — 1.5 |0|. Surprisingly, 
the longitudinal velocity correlations are strongly nega- 
tive for all r and all accessible densities: the velocities 
are anti-correlated. Cm decays slower than exponentially, 
and in the range r = [2a, 16a] can be reasonably well 
described by Cii(r) oc r~ 2 . The structure visible between 
la and 2a is stronger at high density, and its origin is 
not understood. The transverse component Cj_ is much 
smaller and has a shorter range than C | (note the differ- 
ence in vertical scale) . Nearly identical correlations were 
observed without the lid in place. We have performed 
molecular dynamics simulations that closely match the 
conditions of the experiment |12 and find qualitatively 
similar correlations. 

The existence of long-range velocity correlations in a 
randomly forced granular gas has been explained theo- 
retically using a hydrodynamic approach Q. The ve- 
locity correlations result from the balance between the 
excitation of hydrodynamic modes by the random ex- 
ternal forcing and their dissipation through diffusion. 
At shorter range, velocity correlations arise from an in- 
creased probability of re-collision due to the forcing ||. 
This suggests that the correlations could be quite sensi- 
tive to the particular mechanism of energy injection, and 
we believe that this is the origin of the dramatic differ- 
ence between the negative correlations observed in our 
experiment and the positive correlations observed in the 
randomly forced model systems {§ 0, [§ U ■ In those sys- 
tems, energy is injected through white noise forcing of 
each particle independently, while on the smooth plate 
energy is injected only into the vertical motion of the 
spheres and is then transferred to the horizontal motion 
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FIG. 2: (v(0) ■ v(r)) (as described in the text) normalized 
by T for T = 1.5 and p = 0.4 (•), p = 0.5 (o). 



through inter-particle collisions. As a result of this trans- 
fer, collisions will result, on average, in an increase in the 
magnitude of the relative horizontal velocities v r of pairs 
of particles. This effect can be measured directly in the 
molecular dynamics simulations, and at p — 0.5 we find 
that the average value of (v 2 ) for pairs of particles leav- 
ing a collision is 13% higher than for pairs entering a 
collision. Experimental support for this scenario can be 
found in Fig. ^, which shows a calculation of (v(0) • v(r)), 
where the average includes only particles restricted to a 
narrow band along the direction of motion of the par- 
ticle at the origin, either ahead (for positive r) or be- 
hind (for negative r). There is a clear asymmetry with 
respect to zero, showing that balls moving away from 
each other contribute significantly more to the velocity 
anti-correlation than do particles moving towards each 
other. The velocity anti-correlations would presumably 
dissipate through diffusion as the momentum gets trans- 
ferred to surrounding particles through collisions, thus it 
is not surprising that the observed correlations extend 
over several mean free paths. This mechanism for the 
generation of anti-correlated velocities may also operate 
in inhomogeneously forced granular media when the av- 
erage kinetic energy of the grains is anisotropic. 

To test our hypothesis about the origin of the corre- 
lations, we have changed the forcing by using a rough 
surface, so that the vibrating plate also injects energy 
directly into the horizontal motion of the spheres. This 
energy input mechanism is closer to the white noise forc- 
ing of Ref . j7j , and allows for a more direct comparison 
between experiment and theory [Q. 

The velocity correlations measured on the rough plate 
are shown on Fig. ^ for T = 1.5 filjl , and two extremal 
densities. We find that both C\\ and C± are positive 
for all measured p. (v(0) • v(r)), calculated as described 
above, is positive everywhere and is roughly symmetric 
about zero. Molecular dynamics simulations of the layer 
on the rough plate |L2| show similar correlations, and the 
effect of the collisions is to decrease the average value 
of (v 2 .) by 3.7% at p = 0.5. Thus in the absence of 
large vertical to horizontal collisional energy transfer, the 
velocity anti-correlations are replaced by positive velocity 
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Ta(p)f\\ t x(r) (Fig. ||). The functions f\\ and fx have 
different r dependence, but the function a(p) is the same 



for both Cm and C 



FIG. 3: Longitudinal velocity correlations normalized by the 
granular temperature on the rough plate: T = 1.5, and p = 
0.125 (•), p = 0.6 (o). Inset: transverse velocity correlations. 
Data are averaged over 16000 pairs of images (~ 80 balls per 
image) for p = 0.125, and 8000 pairs (~ 420 balls) for p = 0.6. 
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FIG. 4: Scaled longitudinal and transverse velocity correla- 
tions normalized by T. C\\/(aT) (symbols with solid lines) 
and C±/(aT) (symbols with dashed lines) for F — 1.5, 
p = 0.125 (o), p = 0.5 (x), and p = 0.6 (•). Inset: Log-linear 
plot of the scaled velocity correlations (the straight lines are 
exponentials with decay lengths of 2a and 0.6<r.) 



correlations, presumably arising from the random forcing 
as described in Refs. [Q, |sj . The dramatic difference in 
the correlations between the smooth and rough plates 
is a clear indication that the forcing, rather than the 
inelastic interparticlc collisions, determines the velocity 
correlations. 

Interestingly, we find that C\\ and C± can be roughly 
scaled on two different curves of the form Cm j_ = 



/ii and fx are reasonably well de- 
scribed by f(r) oc e~ r / r °, with r a — la for f\\ and 0.6a 
for fx- ft is surprising that this decay length has lit- 
tle density dependence, despite the fact that the mean 
free path estimated from Enskog-Boltzmann kinetic the- 
ory varies from 2.3a for p = 0.125 to 0.16a for p = 0.6. 
A similar scaling relation is predicted in the randomly 
forced model at large r 0, but the scaling is the same 
for both C|| and Cx, and f(r) decays algebraically. 

An important difference between our experiment and 
the randomly forced model is revealed by the velocity dis- 
tributions P(v) of the balls on the rough plate. Figure ^ 
shows the normalized P(v x ) obtained for T — 1.5, and p 
ranging from 0.014 to 0.6. At the lowest density, P(v x ) 
is well described by a gaussian. As p increases, the tails 
of the distribution rise significantly above a gaussian. At 
high density, the velocity distribution resembles those re- 
ported previously on flat plates @, ||. To better under- 
stand the effect of density on P(v), we have followed the 
trajectories of individual balls using a fast camera (838 
frames/s, Q), and measured the time-dependent veloc- 
ity autocorrelation function C v (t) = (v(0) • v(i)) and 
find that C v (t) ~ e _ */ T ^ p - ) at all densities. The scattering 
length defined by I = \/Tt(p) varies from about 0.5er at 
p = 0.6 to about 2.5cr at p = 0.014. At the lowest density, 
I is much smaller than the mean free path for ball-ball 
scattering (~ 20 a) obtained from kinetic theory. This 
suggests that it is ball-plate scattering that is responsible 
for the decay of the velocity auto-correlation function, 
and that the scattering must therefore have a viscous 
component. The evolution of P(v) for a single ball on 
the rough plate can thus be very likely modelled with 
an equation including only random white noise forcing 
and a viscous term, such as the Fokker-Plank equation 
for Brownian motion, which produces a gaussian P(v). 
When the random forcing is instead balanced by ball- 
ball scattering in a modified Boltzmann equation, there 
is an over-population of the tails of the distributions, 
P(v) ~ exp(-|v| 3 / 2 ) J|. This is consistent with the 
experimentally measured P(v) at high densities, where 
the scattering is dominated by ball-ball collisions. At in- 
termediate densities both scattering processes contribute 
and the evolution of P(v) with density can likely be de- 
scribed with a Boltzmann equation consisting of random 
white noise forcing, collisional scattering, and a density 
independent viscous drag. This model has been recently 
investigated |L5|], although primarily at higher inelastici- 
ties, where significant clustering is observed. 

The ball-plate scattering may also explain the lack of 
significant density dependence of the decay length of the 
velocity correlations. We speculate that the density inde- 
pendent scattering with the plate, which unlike ball-ball 
collisions does not conserve momentum, is more effective 
at destroying the long-range velocity correlations, and 
therefore controls the observed decay. Support for this 
scenario is provided by the fact that the decay length for 
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FIG. 5: Normalized velocity distributions P(v x ), on the rough 
plate for T = 1.5 and p = 0.014 (A), p = 0.125 (o), p = 0.6 (•). 
The solid line is a gaussian. (Results for P(v y ) are identical). 



the exponentially decaying Ci | (r) at measured densi- 
ties is close to the ball-plate scattering length deduced 
from the individual trajectories at very low densities (2a 
versus 2.5a). However, a full quantitative understand- 
ing of the velocity correlations observed on the rough 
plate will require a theoretical analysis along the lines of 
Rcf. [Q, including the effect of viscous drag. Our results 
provide a clear experimental demonstration that the ve- 
locity distribution and correlations in a granular gas are 
not determined exclusively by internal dynamics (inelas- 
tic ball-ball collisions) but by a balance between internal 
dynamics and external forcing. 
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